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FOREWORD 


This  Interim  scientific  report  was  prepared  by  Prof.  Jack  L.  Kerrebrock 
of  the  Gas  Turbine  and  Plasma  Dynamics  Laboratory,  Department  of  Aeronautics 
and  Astronautics,  Massachusetts  Institute  of  Technology.  The  work  was 
accomplished  during  the  period  1 March  1976  to  28  February  1977. 

This  work  Is  being  administered  under  the  direction  of  the  Air  Force 
Office  of  Scientific  Research,  supervised  by  Dr.  Bernard  T.  Wolf son.  Program 
Manager . 


I.  INTRODUCTION* 


The  MHO  laser  Is  a nonequilibrium  generator  In  which  the  disequilibrium 
between  the  electron  gas  and  the  translational  degrees  of  freedom  of  the 
carrier  gas  Is  used  to  produce  a population  Inversion  In  a molecular  additive, 
such  as  CO^*  As  the  difference  between  the  electron  and  gas  temperatures  can 
be  as  large  as  3000'’K,  and  the  power  density,  measured  by  the  Joule  heating 
rate,  also  very  large  (up  to  500  watts/cm^  In  generators),  there  Is  the 
possibility  of  gain  and  very  high  lasing  power  densities. 

There  are  several  appealing  features  of  this  concept.  As  the  electric 

power  generated  by  the  MHD  process  Is  delivered  locally  to  the  molecular  gas, 

entirely  within  the  flow,  there  Is  no  need  to  take  high  powers  out  of  the 

power  source,  and  In  turn  put  it  back  Into  a gas  discharge.  This  eliminates 

many  uncertainties  stemming  from  power  conditioning  and  from  electrode 

< 

phenomena.  Further,  while  a high  speed  gas  flow  may  have  to  be  provided  in 
an  electric  discharge  laser  In  order  to  remove  the  energy  rejected  to  the 
lower  laser  level,  this  flow  Is  Intrinsic  to  the  MHD  laser.  The  process  for 
inversion  production  being  essentially  local,  there  Is  a chance  of  achieving 
more  uniform  gain,  hence  better  beam  quality  than  In  externally  excited 
devices.  Consideration  of  the  potential  efficiencies  of  MHD  lasers  must 
recognize  that  It  converts  thermal  energy  Into  optical  energy.  Various 
estimates  have  placed  the  efficiency  of  such  conversion  between  3 and  5 percent. 
This  Is  much  higher  than  Is  achieved  by  gas  dynamic  lasers.  Although  gas 
discharge  lasers  may  approach  20  percent  efficiency  of  conversion  from 

* The  serious  reader  will  note  that  the  Initial  paragraphs  of  this 
introduction  are  Identical  to  those  of  Ref.  1.  It  Is  felt  they 
are  needed  because  of  the  arcane  nature  of  the  subject. 

New  aspects  of  the  work  are  discussed  on  pp.  3-5. 


electric  power  to  optical  power,  this  efficiency  must  be  multiplied  by  that 
of  the  electric  power  system  (and  power  conditioning)  for  which  25  percent 
seems  optimistic.  Thus,  the  overall  efficiencies  of  MHD  and  gas  discharge 
lasers  may  be  comparable. 

The  research  reported  here  is  an  outgrowth  of  research  on  high  power 

density  nonequilibrium  MHD  generators  at  MIT.  This  research  has  led  to  the 

demonstration  that  very  high  levels  of  dissipation  can  be  achieved  In  super- 

-3 

sonic  nonequilibrium  generators,  up  to  500  watts  cm  , with  differences 

2 

between  electron  temperature  and  gas  temperature  approaching  3000”K. 

In  a parallel  program.  It  was  also  demonstrated  that  high  levels  of 

vibrational  excitation,  near  equilibrium  with  the  electron  gas,  can  be 

achieved  In  such  generators.  Excitation  of  carbon  monoxide  was  demonstrated 

3 

In  the  MIT  generator  In  1969. 

The  possibility  of  direct  excitation  of  molecular  gases  to  produce  a 
single-cavity  high  power  density  laser  Immediately  suggested  Itself.  Such  a 
concept,  using  carbon  dioxide  In  a cesium-helium  plasma  has  been  analyzed 
In  Ref.  4. 

The  analysis  and  experience  with  nonequlllbrlum  plasmas  suggested  that 
the  physical  phenomena  of  such  a device  might  be  exceedingly  complex,  since 
the  kinetic  and  thermoehemlcal  complexity  of  molecular  excitation  Is  added 
to  the  already  Intricate  behavior  of  the  electrothermally  unstable  non- 
equilibrium  plasma.  Accordingly,  the  present  study  was  undertaken,  of  the 
behavior  of  a nonequlllbrlum  plasma,  with  mlecular  addition.  In  a strong 
magnetic  field,  but  In  a pulsed  experlswnt  which  stimulates  the  kinetic 
situation  of  an  MHD  laser  without  the  coaq^lexlty  and  expense  of  a high  speed 
flow  system. 
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The  technique  of  modeling  the  flow-induced  electric  field  of  an  MHD 
laser  by  an  applied  electric  field  has  the  great  advantage  that  It  enables 
the  uncoupling  of  fluid  dynamic  and  bulk  plasma  phenomena.  This  technique 
was  used  to  obtain  the  first  experimental  data  on  electrothermal  (Ionization) 
Instabilities  In  nonequilibrium  plasmas,^  and  this  data  has  since  been 
correlated  by  theory^  and  corroborated  by  measurements  In  nonequilibrium 
generators.^  A further  advantage  of  the  pulsed  technique  Is  that  the  very 
complex  kinetic  processes  which  occur  In  the  plasma  as  it  flows  downstream 
In  the  actual  laser,  are  displayed  as  a (temporal)  transient  In  the  pulsed 
experiment,  vastly  simplifying  the  diagnostic  problem. 

Preliminary  experiments  carried  out  by  this  technique  were  reported 
In  Ref.  6.  These  experiments  demonstrated  the  usefulness  of  the  pulsed 
technique,  and  also  demonstrated  gain  and  net  power  output  from  a Cs  seeded 
mixture  of  He,  CO2  and  N2  at  electron  densities  and  current  densities 

appropriate  to  MHD  laser  operation. 

9 

A theoretical  treatment  of  the  physical  processes  In  the  MHD  laser 
produced  results  In  qualitative  agreement  with  the  experiments,  but  over- 
estimated the  Importance  of  deexcltatlon  of  Cs  by  C02>  leading  to  a 
pessimistic  estimate  of  performance. 

The  purposes  of  the  work  described  here  are:  1)  to  carry  out  a systematic 
experimental  study  of  the  gain  and  optical  properties  of  plasmas  suitable  for 
MHD  lasers,  and  2)  based  on  this  Information  to  carry  out  a feasibility 
assessment  of  the  MHD  laser  concept. 

Since  the  last  reporting  of  this  work  In  Ref.  1,  the  results  of  a major 
study  of  MHD  lasers  have  been  reported  In  Ref.  10.  A theoretical  evaluation 
of  the  MHD  laser  concept  Is  presented,  along  with  experimental  results  from 
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flow  systemB  In  which  U x ^ Induced  gain  was  observed  In  Cs  seeded  He-C02 
plasmas.  The  results  were  very  erratic,  evidently  because  of  chemical 


r 

( . 
I . 


^ reaction  between  the  Cs  and  the  C02<  and  difficulty  with  the  salt  windows, 

which  were  chemically  attacked.  This  led  the  Investigators  to  substitute 

Xe  for  Cs  as  a seed;  they  acquired  more  consistent  data  but  gain  only  with 

applied  electric  fields,  not  with  U x B excitation. 

There  Is  general  agreement  that  the  MHD  laser  requires  a plasma 
-2  -5 

composition  of  about  10  CO2  and  10  Cs  In  He  for  successful  operation. 

[The  CO2  fraction  Is  limited  by  the  Joule  heating  available  to  raise  the 

electron  temperature  to  the  levels  required  to  pump  the  upper  laser  level 

I (around  2500”K) . The  Cs  Is  needed  to  provide  the  required  electron  density 

at  these  electron  temperatures.  This  gas  composition  poses  two  problems 

which  are  unique  to  the  MHD  laser  among  CO2  lasers. 

First,  the  low  CO2  concentration  and  absence  of  N2  alter  the  pumping 

and  deactivation  mechanisms.  The  upper  level  pumping  Is  directly  by  electrons 

4 

■ and  is  well  understood.  In  most  treatments  of  CO2  lasers  the  two  lower 

' levels  (100  and  010)  have  generally  been  assumed  to  be  tightly  coupled  and 

j deactivated  by  atomic  collisions  (by  He  In  the  MHD  laser).  We  find  this 

assumption  iu  Invalid  at  the  low  CO2  concentrations  required  by  the  MHD 
laser  because  the  lower  laser  level  (100)  Is  deactivated  by  collision  with 
j C02,  the  rate  being  limiting.  Thus  a three-level  model  Is  required  for 

! the  CO2.  Such  has  been  included  in  the  MHD  laser  model  to  be  discussed  In 

Section  II.  It  Is  shown  there  that  the  (100  010)  transfer  rate  Is  critical 

* to  MHD  laser  performance.  If  lasing  at  10.6  p Is  required.  It  may  be  that 

the  MHD  laser  should  operate  at  9.6  y,  but  this  question  has  not  been 
explored  In  depth  as  yet. 
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Secondly,  the  low  Cs  fraction  and  chemical  reactivity  of  Cs  and  CO^ 
raise  the  question  whether  these  two  species  can  In  fact  coexist  In  the 


plasma.  Early  estimates  of  the  reaction  rate  suggested  that  they  could  for 
the  flow  times  required  In  full  scale  MHD  lasers  such  as  were  projected  In 
Ref.  4,  but  the  reaction  has  led  to  great  uncertainties  In  the  results  of 
all  small  scale  feasibility  experiments.  Indeed  a great  deal  of  the  effort 
of  the  research  described  here  has  gone  to  developing  the  gas  handling 
system  for  the  pulsed  apparatus  so  that  the  chemical  reaction  between  Cs  and 
CO^  can  be  controlled,  and  that  the  amount  of  Cs  present  during  the  experiment 
can  be  accurately  determined.  Section  III  presents  the  methods  evolved, 
which  do  work. 

Section  IV  presents  experimental  results,  which  Include  the  gain  of  the 
MHD  laser  plasma  as  a function  of  current  density  and  magnetic  field,  and 
transient  data  from  which  the  (100-010)  transfer  rate  In  CO2  can  be  Inferred. 
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II.  MHO  LASER  MODELING 


2.1  Summary 

A computer  program  describing  the  performance  of  an  MHO  laser  has  become 
operational.  The  program  combines  a previously  developed  model  of  the  fluid 
flow  In  a nonequilibrium  MHO  generator  with  a new  description  of  CO^  molecular 
kinetics.  The  results  of  this  program  Indicate  that  In  a laser  containing 
one  percent  or  less  of  CO^  the  optical  power  extraction  rate  Is  limited  by 
the  rate  at  which  the  10*0  symmetric  stretch  level  Is  deactivated  to  the  01 '0 
bending  level  by  colllslonal  processes.  This  relaxation  proceeds  an  order  of 
magnitude  more  slowly  than  the  deactivation  of  the  01 '0  level  to  the  ground 
state,  which  had  previously  been  assumed  to  be  the  rate-llmltlng  step. 

2 . 2 Discussion  of  CO^  Kinetics 

Previous  analyses  of  MHD  lasers  employing  CO^  as  the  optically  active 
molecular  species  have  assumed  that  the  symmetric  stretch  (v^)  and  bending 
(V^)  modes  are  In  equilibrium  at  a single  vibrational  temperature  due  to  the 
Fermi  resonance  that  exists  between  them.  As  It  Is  well  known  that  the  rate 
of  deactivation  of  the  01 *0  via  collisions  with  helium  atoms  Is  very  rapid 
(-5.0  X lO^atm  ^sec  ^),  such  a strong  coupling  between  the  modes  Implies  a 
rapid  removal  of  molecules  from  the  10*0  lower  laser  level.  This  has 
Important  consequences  during  the  extraction  of  optical  power  from  the  laser 
gain  medium,  because  It  establishes  an  upper  limit  on  the  rate  at  which 
stimulated  transitions  between  the  00*1  and  10*0  levels  can  occur.  If  this 
rate  of  production  of  10.6  y photons  exceeds  the  rate  of  removal  of  molecules 
from  the  lower  laser  level,  the  population  of  the  lower  level  will  build  up 
very  rapidly  until  the  population  Inversion  Is  destroyed.  Then  the  optical 
power  that  can  be  removed  from  the  laser,  given  by  the  energy  per  transition 
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times  the  photon  production  rate,  is  limited  by  the  colllslonal  deactivation 
rate  of  the  lower  laser  level. 

For  conditions  under  which  the  assumption  of  strong  coupling  between 

symmetric  stretch  and  bending  modes  is  valid,  the  known  rapid  deactivation 

rate  of  the  01 '0  level  leads  to  a prediction  of  high  optical  power  output. 

Thla  prad Ictlon  Is  confirmed  by  the  high  continuous  wave  power  obtained  from 

CO^  - N2  - He  gas  dynamic  and  electric  discharge  lasers.  The  assumption  of 

equilibrium  between  modes  and  has  always  seemed  valid  because  the 

13 

performance  of  these  lasers  closely  matched  prediction,  and  because  of  the 
measurement  of  this  rate  by  Rhodes,  Kelly  and  Javan,  which  Implied  a very 
strong  coupling  (“3.0  x 10®atm  ^sec 

However,  more  recent  studies  have  reported  a variety  of  slower  coupling 
15-22 

rates,  most  of  them  slower  by  an  order  of  magnitude  than  the  RKJ  (Rhodes, 

19 

Kelly,  Javan)  rate.  Bulthuls  performed  experimental  studies  of  the  coupling 

process  In  a CO2  - N2  - H2O  laser  plasma  which  Indicated  that  equilibrium 

between  v.  and  v.  existed  If  the  CO.  concentration  exceeded  12%  of  the  H.O 
1 2 1 2 2 

concentration.  From  this  he  concluded  that  the  primary  coupling  mechanism  was 
C02(10*0)  + C02(00‘’0)  t 2C02(01’0)  + 54  cm‘^ 

20 

These  studies  were  extended  to  CO2  “ N2  - He  plasmas  by  Stark,  by  Murray, 

22  21 
Mltchner,  and  Kruger,  and  by  Jacobs,  Pettipiece,  and  Thome  1.  Based  on 

momentum  transfer  considerations,  Stark  argues  that  the  relative  effectiveness 

of  collision  partners  for  the  10*0  - Ol'O  deactivation  process  is  1:0.46:0.054 


for  CO2  : N2  : He.  Using  this  estimate,  Murray,  Mltchner  and  Kruger  report  a 

value  for  the  CO2  - CO2  rate  of  9.1  x lO^atm  ^sec  while  the  CO2  - He  rate 

Is  4.‘J  X lO^atm  ^sec  If  only  CO2  - CO2  collisions  are  effective  in  this 

18 

process,  as  Rosser,  Hoag  and  Gerry  report,  the  CO2  - CO2  rate  is  slightly 
higher,  1.52  x lO^atm  ^sec 
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Finally,  Jacobs,  Pettlplece  and  Thomas  report  that  the  deactivation 

process  actually  proceeds  via  two  or  more  steps,  the  most  Important  of  which 

are:  ^ 

CO^dO'O)  + M ^ 002(02*0)  + M + 102.8  cm~^ 

002(02*0)  + M ♦ 002(01'0)  + M + 50.1  cm"^ 

where  M Is  either  002(00*0)  or  He.  They  report  a very  rapid  rate  for  the 
8 “1  "1 

first  process  (=<10  atm  sec  for  both  OO2  - OO2  and  OO2  - He  collisions). 

However,  the  second  step  of  the  process  was  found  to  have  a rate  about  ten 

times  slower  (^lO^atm  ^sec  ^).  This  Is  of  the  same  order  of  magnitude  as  the 

rate  determined  by  Seeber^^  for  the  10*0  - Ol'O  relaxation,  and  within  a 

19 

factor  of  ten  of  the  rate  determined  by  Bulthuls  for  this  process  In  a 
CO2  - N2  - H2O  mixture.  It  is  60  to  100  times  larger  than  the  rate  given  by 
Murray,  Mltchner  and  Kruger,  but  10  times  slower  than  that  of  Rhodes,  Kelly 
and  Javan. 

The  first  conclusion  that  may  be  drawn  Is  that  the  detailed  kinetics  of 
the  colllslonal  relaxation  process  of  the  lower  laser  level  C02(10°0)  are 
still  not  well  understood.  There  Is,  however,  a wide  body  of  theoretical 
and  experimental  evidence  to  suggest  that  the  relaxation  rate  Is  from  one  to 
three  orders  of  magnitude  slower  than  previously  believed.  A preliminary 
measurement  from  the  present  work  supports  this  conclusion.  It  and  the  other 

data  are  listed  in  Table  1.  The  evidence  also  Indicates  that  the  CO2-  CO2 
collisions  are  much  more  effective  In  this  relaxation  process  than  either 
CO2  “ He  or  CO2  - H2O  collisions.  Then,  for  low  CO2  concentrations,  the 
relaxation  of  the  01*0  level  by  H2O  or  He  Is  faster  than  the  coupling 
between  the  10*0  and  Ol'O  levels.  This  does  not  pose  any  problems  for  most 
high  power  CO2  lasers,  because  they  generally  operate  with  CO2  mole  fractions 
between  5 and  ISX,  and  H2O/HC  fractions  between  10  and  15X,  the  balance 
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being  made  up  of  N^.  Thus  the  criterion  of  Bulthuls  that  the  CO^  concen- 
tration must  exceed  12Z  of  the  fraction  of  the  01 '0  relaxant  in  order  to  have 
strong  10*0  - 01 *0  coupling  Is  easily  met. 

In  the  MHD  laser  a far  different  situation  prevails.  A study  by 
9 

Lowensteln  Indicates  that  the  CO2  concentration  must  be  kept  below  O.IZ  In 
order  to  avoid  quenching  of  the  excited  electronic  states  of  the  alkali  seed, 
thus  depressing  the  number  density  of  electrons  required  to  excite  CO^ 
molecules  to  the  upper  laser  level.  While  this  low  limit  resulted  from  a 
probable  overestimate  of  the  CO2  - Cs  quenching  cross  section.  Ref.  4 Indicates 
a CO2  fraction  of  about  0.01.  In  this  case,  according  to  Bulthuls'  criterion, 
the  rate-limiting  process  Is  actually  the  10*0  - Ol'O  relaxation  rate. 

All  previous  analyses  of  MHD  lasers  have  assumed  this  rate  to  be  very  fast, 
as  noted  earlier,  so  that  the  populations  of  the  symmetric  stretch  and  bending 
modes  are  In  Boltzmann  equilibrium.  With  this  assumption  the  rate-limiting 
process  Is  naturally  the  relaxation  of  CO2(01'0)  by  He  collisions.  Thus, 
there  Is  the  clear  possibility  that  the  previous  analyses  may  have  over- 
estimated the  optical  power  obtainable  from  an  MHD  laser  by  up  to  an  order 
of  magnitude. 

In  the  present  analysis  the  three  vibrational  modes  of  the  CO2  molecule 
are  treated  Independently,  each  with  its  o%m  vibrational  temperature.  Since 
all  the  other  colllslonal  and  radiative  transfer  rates  between  the  modes  are 
known  to  within  20X,  we  can  study  the  effects  of  the  uncertainty  In  the 
10*0  - Ol'O  relaxation  rate  without  accumulating  too  much  inaccuracy  because 
of  uncertainties  in  the  other  rates.  In  the  present  report  we  illustrate 
the  results  for  the  two  extreM  cases,  namely  the  very  slow  rate  of  Murray, 
Mitchner  and  Kruger,  and  the  very  fast  rate  implied  by  assuming  and  V2 
are  in  equilibrium. 
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the  electron  continuity  (rate)  equation,  the  electron  momentum  equation 

(Ohm's  Lnw),  and  the  electron  energy  equation.  It  is  known  * that  If  the 

Hall  parameter  B exceeds  a critical  value,  the  plasma  becomes  unstable  and 

electrothermal  waves  may  propagate.  These  result  In  temporal  and  spatial 

nonuniformities  In  n and  T which  reduce  the  values  of  6 and  the  conduc- 

e e 

6 26  27 

tlvlty  a.  By  a quasl-llnear  averaging  technique  developed  by  Solbes  * * 

we  may  account  for  these  reductions  by  calculating  effective  values  and 

o^££  for  the  bulk  plasma. 

In  this  quasl-llnear  model  the  electron  number  density  is  first  written 

as  the  sum  of  an  average  term  plus  a small  perturbation: 

n ■ <n  > (1  + n) 
e e 

The  three  electron  equations  are  themselves  averaged  and  combined  to  yield 

an  expression  describing  the  evolution  of  as  the  flow  moves  down  the 

channel.  By  assuming  the  nonuniformities  have  a plane  wave  structure  It  Is 

possible  to  determine  values  of  calculating  the  joule 

dissipation  which  appears  In  the  heat  transfer  term  of  the  electron  energy 

equation.  Manipulation  of  the  averaged  plasma  equation  yields  a second 

2 

differential  equation  describing  the  growth  of  <ri  >,  the  fluctuation  parameter, 
with  flow  distance. 

A previous  numerical  solution  of  the  fluid  equations  for  the  MHD  laser 

11  12 

by  the  United  Aircraft  Research  Laboratory  * assumed  functional  forms  for 

♦ 

a and  6 ...  Furthermore,  constant  values  were  assumed  for  the  critical 
Hall  parameter  (above  which  nonuniformities  develop  In  the  plasma)  and  the 
apparent  Hall  parameter  (which  describes  the  reduction  from  Ideal 

performance  due  to  generator  losing  effects  and  lnterelectroci4^  shorting) . 

In  the  present  analysis  each  of  these  quantities  Is  calculated  explicitly 


at  every  step  In  the  numerical  Integration  using  the  reduction  formulas  due 
6 26  27  24 

to  Solbes  * * and  Cole.  It  Is  felt  that  this  gives  a more  exact 

description  of  the  plasma  state  on  a local  level.  The  values  of 
that  result  give  a more  accurate  calculation  of  the  pumping  rate  of  the  CO^ 
molecules  by  electrons. 

2.3  Coupling  of  MHD  Model  with  CO^  Kinetics 

The  new  kinetic  model  of  the  CO2  molecule  was  combined  with  the  Improved 

MHD  flow  model  to  calculate  the  excitation  rate  of  the  lasing  species  and  the 

optical  power  output.  In  the  present  study  computer  programs  were  developed 

for  both  the  two-temperature  and  three-temperature  CO2  kinetic  isodels,  so  that 

the  effects  of  weak  coupling  between  the  10°0  and  Ol'O  vibrational  levels 

could  be  compared  with  the  equilibrium  result. 

From  a mathematical  viewpoint  the  most  obvious  difference  between  the 

two  models  Is  that  the  three-temperature  version  requires  the  numerical 

Integration  of  seven  ordinary  differential  equations;  the  two-temperature 

model  has  only  six  differential  equations.  In  addition  to  the  two  gasdynamlc 

2 

equations  (for  p^  and  T^)  and  the  two  plasma  equations  (for  <n^>  and  <ri  >), 
each  model  must  solve  a set  of  CO2  vibrational  energy  equations.  For  the 
three-temperature  model  these  are: 


^2e  [ 3 3 3 
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V3-  v^ 


where 


" nv. /k T 
e ^ J -1 


and  represents  the  relaxation  time  of  a molecule  in  mode  1 due  to  a 
collision  with  a member  of  the  j th  species.  In  the  two-temperature  model 
the  first  two  equations  are  combined  as 

i i ” V ■ " [vWh]  * 4;) 


'^3  '^l 


where 


nVj^  2nV2 

®12^V  nv./kT.  ”C0,''‘  nv  /k  T "co, 
e^  ^-1  ^ e ^ ^-1  ^ 

The  equation  for  the  upper  level  energy  remains  unchanged. 

We  now  present  some  results  of  calculations  performed  for  identical 

geometries  and  gas  compositions  using  the  two  different  models.  The  flow  Is 

assumed  first  to  be  expanded  through  a wedge  nozzle  5 cm  In  length  with  a 

divergence  half  angle  of  10.7*  and  an  area  ratio  at  the  exit  ^ 5.6, 

yielding  an  Isentroplc  Mach  number  of  4.  In  this  region  the  fluid  Is 

18  “3 

assumed  to  be  prelonlzed  so  that  n^-  8.5  x 10  m at  the  entrance  to  the 
channel.  The  gas  next  enters  the  wedge-shaped  channel  section  with  a width 
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of  54  cm  and  divergence  half  angle  of  0.3”.  A magnetic  field  of  4.0  tesla  Is 
Initiated  1 cm  Into  the  channel.  Mirrors  of  reflectivity  0.98  and  0.95  are 
positioned  parallel  to  the  channel  axis,  creating  an  optical  axis  perpendic- 
ular to  the  flow  and  magnetic  field. 

The  logic  of  the  program  Is  structured  so  that  Initially  the  optical 
Intensity  is  zero.  At  each  point  In  the  Integration  the  small  signal  gain 
Is  calculated  and  compared  to  the  threshold  gain  for  the  laser  cavity,  given 

‘’y  g - - «.n 

where  ” Mirror  Reflectivities 

L >•  Mirror  Separation  Distance 

When  this  quantity  Is  exceeded  for  the  first  time,  the  calculation  of  the 
Intensity  is  Initiated,  based  on  the  condition  ■ 0,  where 


g * g (I.  T^»  T,  Tj^qq. 

Is  the  local  value  of  the  saturated  gain.  In  this  way  the  extraction  of  power 


from  the  laser  cavity  begins  as  soon  as  the  pimping  rate  of  the  upper  laser 
level  reaches  a value  that  pushes  the  gain  above  the  threshold. 

The  gas  composition  considered  here  has  the  following  mole  fractions: 

«c.  - 

x„.  • 0.98999 


Stagnation  conditions  are  given  by: 

Pq  • 20  atm 
Tq  - 2090”K 


These  conditions  result  In  a mass  flow  rate  of  1.56  kg/sec. 
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Temperature  response  profiles  for  the  two-  and  three-temperature  models 
are  presented  In  Figs.  2.1  and  2.2  respectively.  The  numerical  results  from 
the  calculations  are  given  In  Tables  2 and  3.  Finally,  plots  of  the 
local  optical  power  density  gl  are  Included  as  Fig.  2.3,  giving  a comparison 
between  the  two  models. 

Calculations  with  the  three-temperature  model  use  the  lower  level 
deactivation  rate  of  Murray,  Mltchner  and  Kruger,  which  assumes  that  both 
CO^  and  He  are  Involved  In  the  relaxation  process,  with  rates  given  by: 
’‘lOO-CO^  " lO^tm'^sec’^ 

'^lOO-He  - lO^atm'^sec’^ 

This  set  of  parameters  yields  a somewhat  higher  relaxation  rate  than  that 
given  by  Bulthuls,  who  assumed  that  only  CO2-CO2  collisions  were  effective 
In  the  deactivation  of  the  lower  laser  level.  This  does  not  seem  entirely 
realistic,  although  It  does  Indicate  the  much  greater  effectiveness  of  the 
CO2  molecule  In  the  colllslonal  relaxation  process.  For  the  conditions 
assumed  In  this  simulation,  the  two-temperature  model  gives  a relaxation 
time  for  the  lower  laser  level  of  about  1 microsecond;  the  three-temperature 
model  gives  about  10  microseconds. 

The  Immediate  consequences  of  this  can  be  seen  by  examining  the  curves 
for  Tj^qq  In  Figs.  2.1  and  2.2.  The  vibrational  temperature  (and  hence  the 
population)  of  the  lower  laser  level  remains  hung  up  at  a high  value  during 
and  after  expansion  through  the  nozzle  In  the  three-temperature  model. 

In  fact,  with  the  low  CO2  concentration  employed  In  the  MHD  laser,  the  three- 
temperature  model  Indicates  that  the  lower  laser  level  actually  relaxes  more 
slowly  than  the  upper  (00*1)  level.  A population  Inversion  Is  not  possible 
until  CO2  molecules  get  excited  by  electron  Impact  In  the  MHD  channel  section. 
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Teaperature  response  profiles  for  the  two-  and  three-teaperature  andels 
are  presented  in  Figs.  2.1  and  2.2  respectively.  The  numerical  results  from 
the  calculations  are  given  in  Tables  2 and  3.  Finally,  plots  of  the 
local  optical  power  density  gl  are  Included  as  Fig.  2.3,  giving  a comparison 
between  the  two  models. 

Calculations  with  the  three-teaperature  model  use  the  lower  level 
deactivation  rate  of  Murray,  Mltchner  and  Kruger,  which  assuaes  that  both 
CO^  and  He  are  involved  in  the  relaxation  process,  with  rates  given  by: 
‘'lOO-COj  “ 10^«ta’^sec‘^ 


"lOO-He 


4 -1  -1 

■ 4.8  X 10  ata  sec 


This  set  of  parameters  yields  a somewhat  higher  relaxation  rate  than  that 
given  by  Bulthuis,  who  assumed  that  only  COj-COj  collisions  were  effective 
in  the  deactivation  of  the  lower  laser  level.  This  does  not  seem  entirely 
realistic,  although  it  does  indicate  the  much  greater  effectiveness  of  the 
CO^  molecule  in  the  collisional  relaxation  process.  For  the  conditions 
assumed  in  this  slaulation,  the  two-temperature  model  gives  a relaxation 
time  for  the  lower  laser  level  of  about  1 microsecond;  the  three-teaperature 
model  gives  about  10  microseconds. 

The  immediate  consequences  of  this  can  be  seen  by  examining  the  curves 
for  In  Figs.  2.1  and  2.2.  The  vibrational  teaperature  (and  hence  the 

population)  of  the  lower  laser  level  remains  hung  up  at  a high  value  during 
and  after  expansion  through  the  nossle  in  the  three-teaperature  model. 

In  fact,  with  the  low  CO^  concentration  employed  in  the  MHD  laser,  the  three- 
temperature  model  indicates  that  the  lower  laser  level  actually  relaxes  more 
slowly  than  the  upper  (00*1)  level.  A population  inversion  is  not  possible 
until  CO2  molecules  get  excited  by  electron  impact  in  the  MHD  channel  section. 
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On  the  other  hand.  Table  2 shows  that  the  two-temperature  model 

predicts  some  gain  In  the  nozxle  region  due  to  the  selective  vibrational 

freezing  of  molecules  In  the  upper  laser  level,  the  same  principle  by  which 

a gasdynamlc  laser  operates.  The  low  gain  prediction  of  the  three-temperature 

model,  valid  at  CO^  concentrations  on  the  order  of  one  percent  or  less, 

has  important  Implications  for  high-pressure  CO2  gasdynamlc  lasers  which 

28 

operate  In  a so-called  "C02'atarved"  mode.  Indeed,  measurement  of  small- 
signal  gain  In  a gasdynamlc  laser  with  99Z  He  - IZ  CO2  composition  represents 
one  method  of  checking  the  slow  lower  level  relaxation  rate. 

Since  the  Induced  transition  rate  In  a laser  oscillator  Is  limited  by 
the  deactivation  rate  of  the  lower  laser  level,  we  can  expect  that  with  a 
deactivation  time  of  almost  10  microseconds  the  three- temperature  model  will 
predict  optical  power  production  rates  about  ten  times  less  than  those  given 
by  the  two-temperature  calculations.  This  Is  seen  to  be  the  case  by 
examination  of  Tables  2 and  3,  and  Fig.  2.3.  The  abrupt  bumps  In 
Figs.  2.1  and  2.3  appear  to  have  no  physical  cause;  they  are  believed  to  be 
the  result  of  a numerical  Instability  or  computer  roundoff  error. 

2 . 4 Directions  for  Further  Research 

With  the  computer  simulation  of  HHD  laser  performance  now  operational, 
a program  can  be  carried  out  to  determine  more  optimal  design  features  for 
such  a device.  Thus  far  our  choice  of  operating  conditions  has  been  rather 
arbitrary.  We  are  about  to  undertake  a parametric  study  of  the  effects  of 
channel  geometry,  CO2  concentration,  magnetic  field,  and  mirror  reflectivities 
on  the  perTormance  of  an  MHD  laser. 

However,  It  appears  that  the  two  most  critical  factors  are  the  CO2  (10*0) 
relaxation  rate  and  the  cross  section  for  quenching  of  the  6Pj^y2  "tate  of  Cs 
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by  CO^  encounters.  Each  of  these  quantities  will  be  treated  as  a parameter  In 
the  computer  program.  Several  runs  will  be  made  over  a range  of  values  of 
each  quantity  to  determine  which.  If  any,  permit  the  practical  operation  of 
an  MHD  laser. 

It  would  clearly  be  desirable  to  find  a gas  additive  that  does  not 
seriously  upset  the  energy  balance  In  the  non-equlllbrlum  MHD  plasma,  and 
which  rapidly  relaxes  the  CO2  lower  laser  level.  For  example.  It  has  been 
reported^®  that  the  deactivation  rate  of  CO2  (10*0)  by  Xe  is  of  the  same  order 
of  magnitude  as  the  rate  due  to  CO2  Itself. 
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III.  EXPERIMENTAL  PROCEDURE 
3.1  The  Pulsed  Simulation  Technique 

The  pulsed  experiment  simulates  by  a transient  the  convected  history  of 
the  fluid  as  it  passes  through  the  high  speed  flow  of  an  MHD  laser.  In  earlier 
phases  of  the  present  work,  the  transient  was  used  only  In  the  electrical 
conditions.  A steady  flow  of  He  + Cs  + CO2  was  established  through  a quartz 
test  section.  A pulsed  electric  field  and  magnetic  field  then  established  the 
MHD  plusm.'!  conditions  desired.  Experience  with  this  system  showed  that  the 
Cs  concentration  could  not  be  controlled  well  enough  to  do  a systematic 
experiment.  The  reaction  of  Cs  and  CO^  made  the  actual  Cs  concentration  very 
uncertain.  Kurthermore,  the  continuous  flow  of  He  + Cs  -t*  CO2  through  the  test 
section  quickly  coated  the  diagnostic  windows  and  mirrors,  requiring  frequent 
disassembly  for  cleaning. 

Accordingly  the  gas  handling  system  was  redesigned  for  pulsed  operation. 
The  sequence  of  events  which  occur  during  a test  Is  shown  schematically  In 
Pig.  3.1.  Helium  Is  stored  In  a heated  flask  containing  copper  wool  wet  with  Cs. 
VHien  a valve  is  opened  this  mixture  blows  down  to  vacuum  through  the  test 
chamber  giving  the  pressure  transient  sketched.  Near  the  peak  pressure, 
the  magnetic  field  is  pulsed  on,  and  the  shutter  on  the  CO2  probe  laser  Is 
opened.  Finally,  the  electric  field  Is  pulsed  to  pass  a prescribed  current 
through  the  plasma  after  a short  electrostatic  prelonlzer  pulse. 

The  critical  diagnostics  consist  of: 

1.  The  CO2  probe  laser  and  detector. 

2.  The  Cs  absorption  measurement. 

3.  Electron  density  measurement  by  the  Cs  bound-free  continuum. 
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A.  CO  measurement  by  line  absorption. 

5.  Hall  and  axial  fields  by  probes  In  contact  with  the  plasma. 

Some  unusual  features  of  the  apparatus  will  be  discussed  In  the  following. 

3 . 2 Pulsed  Flow  System 

The  riow  system  Is  sketched  In  Pig.  3.2.  Before  a series  of  tests, 
liquid  Cs  is  flushed  through  the  right  hand  stainless  steel  vessel  (of  about 
0.3  9,  capacity)  so  as  to  thoroughly  wet  the  Cu  wool.  The  liquid  Cs  Is  then 
returned  to  the  left  vessel.  During  operation  the  pneumatic  valve  Is  opened 
allowing  the  He  saturated  with  Cs  to  flow  through  the  test  section  for  a 
period  of  about  700  ms,  the  maximum  pressure  being  about  63  mmHg.  This  provides 
reliable  seeding  with  Cs.  It  also  minimizes  fogging  of  the  test  section 
windows. 

To  minimize  the  Cs-CO^  chemical  reaction,  the  CO^  Is  Introduced  as  shown,  ^ 

Into  the  flowing  He-Cs  mixture,  the  operation  of  the  solenoid  valve  being  j 

timed  to  coincide  with  that  of  the  pneumatic  valve  so  as  to  give  a nearly  | 

constant  CO^  fraction. 

3. 3 Cs  Concentration  Measurement 

Determination  of  the  actual  concentration  of  metallic  Cs  in  the  discharge 
has  been  especially  troublesome  in  MHD  laser  investigations.  Here  the  measure- 

O 

mcnt  Is  done  by  resonance  absorption  at  4550  A of  light  from  a Cs  lamp. 

To  minimize  noise  and  record  a reference  on  each  scope  trace  the  light  from 
the  lamp  is  chopped  at  about  15  kHz,  giving  ouiny  periods  during  a gas  flow 
pulse,  as  Indicated  In  Pig.  3.1.  The  change  In  height  of  the  pulses  gives 
directly  Che  absorption  by  the  Cs  In  the  tube,  and  since  the  record  spans  the 
whole  gas  flow  time,  any  fogging  of  the  windows  can  be  detected. 

The  Cs  absorption  measurement  was  calibrated  by  means  of  a quartz  test 
chanber  scaled  off  with  60  mmHg  of  He  and  heated  In  an  oven. 
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Tn  the  experiments  to  be  described  below,  the  discharge  occurred  late  in 


the  gas  flow  pulse,  as  shown  In  Fig.  3.1,  where  the  Cs  absorption  was  very 
small,  so  that  the  Cs  concentration  was  not  directly  measured.  Estimates  of 
Its  value  were  Inferred  from  other  data,  as  will  be  explained. 

3.4  CO^  Probe  Laser 

The  measurement  of  gain  Is  done  by  means  of  a three-pass  probe  laser 
beam  at  10.6  p as  shown  in  Fig.  3.3.  A shutter  Is  opened  just  before  the 
electric  field  pulse  to  establish  the  detector  zero  for  each  test.  As  the 
detector  Is  AC  coupled,  the  gain  signal  Is  superimposed  on  a decay  as  sketched 
in  Fig.  3.1.  The  overall  time  resolution  of  the  system  Is  better  than  0.1  ps, 
a point  which  will  be  Important  In  Interpretation  of  the  gain  data. 
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IV.  EXPERIMENTAL  RESULTS 

With  the  above  described  modifications  to  the  appratus,  consistent  gain 

measurements  have  been  achieved  over  a range  of  conditions  appropriate  to  MHD 

lasing  operation.  The  results  are  conveniently  divided  Into  those  '>btained 

without,  and  with  magnetic  field. 

A . 1 Mcusuremcnts  Without  Magnetic  Field 

To  chock  the  flow  system,  discharge  circuits  and  Instrumentation,  a series 

of  experiments  was  first  conducted  without  CO^t  l.e.  with  a Cs  seeded  He  plasma. 

-2  -1 

The  current  was  about  0.5  amp  cm  and  the  conductivity  0.093  mho  cm  . The  Cs 

-4 

absorption  measurement  Indicated  a Cs  mola  fraction  of  1.0  x 10  at  the 
pressure  peak.  If  there  were  no  mechanisms  for  removal  of  Cs  from  the  mixture, 
we  would  expect  this  same  mole  fraction  of  Cs  In  the  plasma  at  the  time  of  the 
discharge  some  100  ms  later.  Analysis  of  the  conductivity  results  Indicated 
that  the  Cs  mole  fraction  was  In  fact  closer  to  2 x 10  This  calculation 
consisted  essentially  In  estimating  the  electron  temperature,  then  inferring 
the  Cs  concentration  from  the  conductivity,  which  depends  principally  on  the 
electron  density. 

There  is  therefore  some  mechanism  acting  to  remove  the  metallic  Cs  from 
the  plasma,  the  time  constant  for  removal  being  of  the  order  of  35  ms. 

Whether  this  mechanism  Is  oxidation  by  traces  of  air,  diffusion  to  the  walls, 
or  some  other  Is  not  known  at  present.  In  any  case.  It  seems  to  be  reproducible, 
hence  can  be  dealt  with . 

Upon  addition  of  CO^  to  a mole  fraction  of  0.0074  In  the  He-Cs  mixture 
with  a total  pressure  of  53.8  mmHg  , gain  was  measured  as  shown  In  Fig.  3.3 
UH  a function  of  time.  The  principal  features  of  this  time  history  of  the 

t 

gain  are  as  follow: 
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;i.  The  gain  rises  In  a tine  of  about  10  )is  to  a value  of  0.08  cni~^, 

1h  nearly  constant  for  some  100  ys,  then  decreases,  becoming  strongly  negative 
before  the  current  is  terminated  at  330  ys. 

b.  The  negative  gain  (absorption)  then  relaxes  on  about  500  ys  time 

scale. 

It  should  be  noted  first  that  the  current  Is  actually  constant  from  about 

1 ys  on,  the  Initial  rounding  being  due  to  the  measurement  system. 

The  decrease  of  gain  from  100  ys  on  Is  believed  to  be  due  to  pumping  of 

the  lower  laser  level  (100)  by  electrons.  As  this  level  Is  heavily  populated. 

It  absorbs  the  10.6  y probe  beam.  When  the  current  Is  terminated,  pumping  of 

this  level  by  electrons  stops  also  and  its  population  decays,  giving  a 

corresponding  reduction  in  absorption. 

This  decay  therefore  yields  a measurement  of  the  deactivation  of  the  100 

level  by  both  100-020  collisions  and  by  the  100-He  mechanism.  The  decay  time 

-4 

being  “50  ys  and  the  CO^  pressure  =5  x 10  atm  we  find  a rate  constant  of 

4 X lO^atm  ^s  which  Is  to  be  compared  to  the  values  given  in  Table  1. 

It  falls  in  the  general  range  of  the  published  values. 

It  should  be  noted  that  measurement  of  the  CO  concentration  by  emission 

at  a line  showed  no  significant  amount  during  the  discharge.  Indicating  that 

CO^  dissociation  did  not  contribute  to  the  gain  termination. 

A series  of  measurements  at  different  current  densities  yielded  peak 

-2 

gains  as  shown  In  Fig.  3.4,  there  being  a maximum  at  about  0.4  amp  cm  for 
this  set  of  gas  conditions.  The  peak  gain  Is  close  to  the  computed  values 
cited  in  Section  TI,  however  larger  values  have  been  measured,  as  noted  below. 
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4.2  Measurements  with  Magnetic  Field 


A set  of  three  experiments  v'lth  different  magnetic  field  strengths  Is 
presented  in  Fig.  3.5,  the  current  and  electric  field  in  the  upper  frames  and 
the  10.6  )j  gain  at  the  bottom.  As  noted  previously,  the  current  Is  actually 
a square  pulse  on  this  time  scale. 

The  axial  electric  field  (in  current  direction)  Is  initially  high  as  the 
CO^  is  pumped,  and  the  gain  rises  on  the  same  time  scale  that  the  electric 
field  falls. 

For  zero  magnetic  field  for  these  conditions  the  gain  then  goes  negative 
on  a 100  MS  time  scale,  then  relaxes  as  in  Fig.  3.3. 

When  the  m.agnctlc  field  Is  Increased  from  zero,  the  gain  displays  a pulsed 
behavior,  the  time  scale  for  the  pulses  being  of  the  order  of  10  ys,  and  the 
amplitude  Increasing  with  magnetic  field.  This  appears  to  be  due  to  the 
Lorentz  force,  which  generates  a pressure  wave  propagating  perpendicular  to 
the  current  and  magnetic  field  and  (probably)  generating  a nonuniform  discharge 
which  carries  the  excitation  out  of  the  path  of  the  gain  probe  beam.  Taking 
the  sound  speed  in  the  tube  as  1000  ms  ^ we  find  a full  wave  period  in  the  2 cm 
cube  of  20  MS,  which  Is  close  to  the  observed  period.  The  ratio,  JxB/p  P.  0.3  m 

Tt  Is  concluded  that  the  peak  values  from  traces  such  as  these  on  Fig.  3.5 
nru  the  gains  which  would  be  realizable  in  a flow  system.  A series  of  such 
peak  values  arc  tabulated  In  Table  4 along  with  other  parameters  of  the 
plasma.  The  peak  gain  values  arc  plotted  in  Fig.  3.6. 

It  will  be  seen  from  Fig.  3.6  that  the  gain  peaks  with  Increasing 
magnetic  field  for  fixed  current  density,  CO2  fraction  and  pressure.  This  is 
probably  due  to  the  pressure  wave  phenomenon  discussed  above.  Since  for  fixed 
current  the  dissipation  and  hence  the  electron  temperature  increase  with 


increasing  magnetic  field,  we  would  expect  the  gain  also  to  increase,  and 
gradually  saturate. 


The  measured  gains  are  quite  large  (by  MHD  laser  standards)  and  in  fact 
larger  than  the  small  signal  gains  cited  in  Section  II.  Reference  to  Tabled 
shows  that  other  plasma  parameters  are  as  expected  for  an  MHD  laser  plasma. 

The  apparent  Hall  parameter  decreases  slightly  with  increasing 

magnetic  field  as  expected  for  an  eletrothermally  unstable  plasma  . 

The  conductivities  and  electron  densities  are  in  the  expected  range,  much 
larger  than  for  usual  gas  discharge  laser  plasmas. 

A critical  question  is  whether  the  electric  field  is  in  the  range 
accessible  by  UxB  action.  For  the  magnetic  field  used  in  these  experiments,  and 
and  a flow  velocity  of  3000  ms  ^ (Mach  number  of  3)  we  estimate 
UxB  = 1500  volt  m ^ * 15  volt  cm  ^ 

which  is  slightly  below  the  measured  values.  Thus  the  magnetic  field  in  an 
MHD  laser  using  the  experimental  plasma  conditions  would  have  to  be  at  least 
1 tesla.  At  state  of  the  art  fields  of  4 to  5 tesla  the  UxB  field  would  be 
several  times  that  required. 
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V.  SUMMARY  AND  CONCLUSIONS 

The  rundamenCal  mechanisms  of  lasing  phenomena  In  MHD  lasers  are  being 
studied  by  a combination  of  numerical  modeling  and  a pulsed  experimental 
simulation  of  the  MHO  flow  laser.  In  the  kinetic  representation  of  the  CO^ 
a three  level  model  Is  used  which  accounts  for  a finite  coupling  between  the 
symmetric  stretch  and  bending  modes.  The  calculations  using  published 
coupling  rates  between  these  (100  -*■  020,010)  levttls  Indicate  that  the  coupling 
may  limit  MHD  laser  output  at  10.6  |i.  It  Is  conjectured  that  this  limit  can 
be  avoided  by  oscillation  at  9.6  p.  The  pulsed  experiment  has  demonstrated 
gain  as  high  as  0.3  percent  per  cm  at  MHD  laser  conditions,  with  gas 
composition  He  + 0.0074  CO2  + 10  ^ Cs  and  total  pressure  50  mmHg.  The  current 

I _2  -1 

density  was  approximately  0.3  amp  cm  , the  electric  field  20  volts  cm  and 
the  magnetic  field  0.5  tesla.  It  Is  concluded  that  plasma  conditions 
> suitable  for  MHD  laser  operation  have  been  achieved  In  the  pulsed  simulation. 

The  numerical  modeling  Indicates  that  these  conditions  can  be  reached  in  an 
MHD  flow  laser. 
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TABLE  1 

SUMMARY  OF  REPORTED  VALUES  OF 
C02(l0“o)-C02(01'0)  RELAXATION  RATE 
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TEMPERATURE  RESPONSE  PROFILES  FOR  THREE-TEMPERATURE 
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